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, Abstract. A complete description of W and Z boson production at high- 

energy colliders requires the resummation of large Sudakov logarithms which 
dominate the production at small transverse momentum. Currently there are 
^ ' two techniques for performing this resummation: impact parameter space and 

' transverse momentum space. We argue that the latter can be formulated in a 

i way which retains the advantages of the former, while at the same time allowing 

, a smooth transition to finite order dominance at high transverse momentum. 



' 1. Introduction 

o' 

The search for a correct description of initial gluon radiation has recently become a 
very active field of theoretical research, stimulated by the large amount of experimental 
data now available. In particular, gluon radiation plays an important role in vector 
boson (W, Z) production at hadron colliders, i.e. processes that provide precision 
■ measurements of Standard Model parameters and appear to be significant backgrounds 

to new physics phenomena. For example, one of the methods of determining Mw 
requires a very accurate description of the W's transverse momentum distribution. At 
small transverse momentum (qr), this distribution is dominated by large logarithms 
\n{Q 2 /q 2 ^) which are directly related to the initial emission of soft gluons. Therefore at 
sufficiently small qx fixed-order perturbation theory breaks down and the logarithms 
must be resummed. 

The origin of the large logarithms is visible already at the leading-order (LO) 
level — the contribution from real emission diagrams for qq — > Vg contains a term of 

the form ^£ In ( 

In the case of higher-order contributions, where more gluons are emitted, 
the logarithmic divergence becomes even stronger. It can be shown that in the 
approximation of soft and collinear gluons with strongly ordered transverse momenta 
kr, i.e. 

k 2 T%i < fc|. 2 < ... < k T . N <q\<^Q 2 (1) 
the dominant contributions to the cross section are 
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where = 2Cp. Due to the hi 2N ~ x (Q 2 / q 2 ^) structure, this approximation is 
commonly known as the Double Leading Logarithm Approximation (DLLA) . Although 
the approximation of truncating the series in (|2|) at finite order clearly breaks down 
once a s hi 2 (Q 2 /q T ) ~ f , the form of the coefficients allows (||) to be resummed, giving 
a so-called Sudakov factor |Q: 




Resummation of (|2|) to all orders gives a finite but unphysically suppressed result 
in the small qx limit. This suppression is caused by the vanishing of strongly-ordered 
phase-space, in which overall transverse momentum conservation is ignored. The 
result in (||) corresponds to a configuration in which a single soft gluon balances the 
W transverse momentum, giving the overall ln(Q 2 / 'q T ) / ' q\ term, while all other gluons 
have transverse momenta <C qr- This is not the dominant configuration in the small 
qx limit. Equally important are non-strongly-ordered contributions corresponding 
to the emission of soft (~ qx) gluons whose transverse momenta add vectorially to 
give the overall q-j- of the W. Although such contributions are formally subleading 
order- by-order, they do dominate the cross section in the region where the Sudakov 
form factor suppresses the (formally) leading DLLA contributions. The non-leading 
'kinematical' logarithms are correctly taken into account by imposing transverse 
momentum conservation (rather than strong ordering) , and this is most easily achieved 
in Fourier transform space (b— space). 



2. b space method 



The b— space method is based on introducing a 2— dimensional impact parameter 
vector b, the Fourier conjugate of qx 0. First, the cross-section (with a delta 
function for transverse momentum conservation included) is Fourier transformed, 
then the conjugated cross section a(b) is computed, and finally this is transformed 
back to momentum space. The advantage of the b— space method is that the soft 
gluon factorisation property is retained even in the presence of transverse momentum 
conservation: 
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This allows for the development of a general expression resumming all terms of the 
perturbation series which are at least as singular as l/q T when qx — > Jl (at the 
parton level): 
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with ao = Ana 2 / (9s), and where 

S(b,Q 2 ) = - C % 



ln[^)A(a s (^)) + B(as(fl 2 )) 
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with bo = 2exp(— 7b). The first two coefficients in each series (0) can be obtained Q 
from the exact fixed-order perturbative calculation in the high qx region by comparing 
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the logarithmic terms therein with the corresponding logarithms generated by the first 
three terms of the expansion of exp(5(6, Q 2 )) in (JsJ) . 

Although the b— space method succeeds in recovering a finite, positive result in the 
qr — ► limit, it suffers from several drawbacks. The first is the difficulty of matching 
the resummed and fixed-order predictions. Since the resummation is performed in the 
Fourier-conjugated space one loses control over which logarithmic terms (in qr — space) 
are taken into account. Therefore there is no unambiguous prescription for matching; 
existing prescriptions require 'unsmooth' switching from resummed to fixed-order 
calculation at some value of qr ■ Secondly, since the integration in (||) extends from 
to oo, it is impossible to make predictions for any qx without having a prescription 
for how to deal with the non-perturbative regime of large b. One prescription is to 
artificially prevent b from reaching large values by replacing it with a new variable 6* 
and by parametrising the non-perturbative large-6 region in terms of a form factor. 
The cross section then reads 

y^e^ / dx A dx B 5{xAXB ~%-)\ [ dbb,J Q (q T b) 
Jo & z Jo 
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f q /A{xA,y) fq/B(xB,^-)F^ b p (Q,b,XA,x B )exp[S(b t ,,Q)} + (q «-> q) 



where, for example, the 'freezing' of b at 6* is achieved by 

b* = , , 6* < b Vlm , 

y/l + itybum) 2 

with the parameter 6i; m ~ 1/Aqcd separating perturbative and non-perturbative 
physics. The detailed form of the non-perturbative function FK remains a matter 
of theoretical dispute (for a review see jq|). In a very simple model in which the non- 
perturbative contribution arises from a Gaussian 'intrinsic' kr distribution, one would 
have F ~ exp(— k6 2 ). The data are not inconsistent with such a form, but suggest 
that the parameter re may have some dependence on Q and x. 



3. qr— space method 

The difficulties mentioned above could in principle be overcome if one had a method of 
performing the calculations directly in transverse momentum space. Given an insight 
into which logarithmic terms are resummed, it should be fairly straightforward to 
perform matching with the fixed-order result. Moreover, the non-perturbative input 
would be required in (and would affect) only the small qr region. 

Three techniques have been proposed for carrying out resummation in 
qr— space || S|. For a detailed discussion of the differences between them the 
reader is referred to J^] . The starting point for all techniques is the general expression 
in impact parameter (b) space for the vector boson transverse momentum distribution 
in the Drell-Yan process ||] at the quark level. In the approach of S , the final result 
(in the simplest case with fixed coupling as and retaining only the leading coefficient 

in the A and B series of (^)) is of the form 

1_ da_ _ ^ ^ L 2 ~ (-2A)( jv - 1 ) ^ / N - 1 
aodql'qf 2 2- (JV-1)! £^ { m 

(8) 
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Here L — \n(Q 2 /q 2 ^), A = a s Cp/T:, and the numbers r TO are denned by [j] 

) 

dy My) ln m if). 
The r m can be calculated explicitly using the generating function 
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(10) 



so that e.g. t = 1, t\ = t 2 — 0, t 3 = — 1£(3) etc. 

Notice that by setting all r m coefficients (except tq) to zero one would immediately 
recover the DLLA form (^). Since there are no explicit subleading logarithms in (||), 
other than those related to kinematics, the presence of the r m coefficients must 
correspond to relaxing the strong-ordering condition. This can be checked explicitly 
by performing the 'exact' 0{a 1 s ) calculation in transverse momentum space: 
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Strong ordering is equivalent to replacing the 8— function by 8^ 2 \kTi ~ <7t)x 
9{kj, 1 — k\ 2 ) +(1 <-> 2). This gives only the leading L 3 term on the right-hand side. 
The C(3) term represents the first appearance of the (kinematic) T3 coefficient of (||). 

In principle the formalism presented above allows for an inclusion of any number 
of such subleading kinematic logarithms, defined by the cut-off value iV max . Figure [l] 
shows that for small values of qt the approximation of the b— space result improves 
with increasing iV max . Therefore by retaining sufficiently many terms one can obtain 
a good approximation (i.e. adequate for phcnomcnological purposes) to the 6— space 
result by summing logarithms directly in qx space. N The technique developed so far 
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Figure 1. The 6-snace result (partem level, fixed coupling, only A' 1 )) compared 
to the expression (p|), calculated for various values of iV max . Here r) = q^/Q 2 ■ 



II The T m are called b m (oo) in 
K However, due to the lack of knowledge of A' 3 ', S' 3 ', etc. it is only possible 
where not more than the first four 'towers' of logarithms are fully resummed, see 



obtain a result 
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can be extended to include subleading A and B coefficients, the running coupling and 
parton distributions, thus yielding a 'realistic' expression for the hadron-level cross 
section. The result is too lengthy to reproduce here, but can be found in ^ 
A comparison between the (Tevatron) DO data ml\ on the Z transverse momentum 
distribution and two of the qx~ space predictions |3, || is shown in Figure |. No non- 
perturbative contribution is included, and in fact the fits can be improved slightly by 
including a modest amount of Gaussian intrinsic kr smearing, as discussed above. 
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Figure 2. Comparison of the theoretical predictions derived in the Ellis- Veseli 
(EV) D and Kulesza-Stirling (KS) g, [l(j approach with DO data Jll] . 



4. Conclusions 

We have discussed two methods of performing resummation in vector boson production 
at hadron colliders. We have argued in favour of the qx— space method, which can be 
extended to systematically take into account subleading logarithms of kinematic origin. 
However more work is required. Although the qr~ space method provides a simple 
matching prescription, the form of the non-perturbative function in this approach 
remains an open theoretical issue. In particular, the current lack of understanding 
of the x and Q dependence of the non-perturbative contribution is a serious limiting 
factor in predicting the q? — > behaviour of the distribution at the LHC. 
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